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We measure the interface tension near the phase transition endpoint 
j^ , of the 3d St/(2)-Higgs model. The tunnel correlation length method is 

used and compared to other approaches. A modified scaling behaviour for 
O ,' the mass gap as function of the transverse area is proposed. 
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^ ■ 1 Introduction 

It is now almost established that the symmetry restoring electroweak phase 
transition changes into a smooth crossover if the mass of the so far elusive Higgs 
particle would be equal to the W-mass (or less within 10 per cent) |l|, g, |[ [|. 
This emerges from lattice studies in the effective 3-dimensional model as well 
as from 4-dimensional simulations, the latter so far only with a rough lattice 
spacing of order 0(1/(2T)) with temperature T. 
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In order to quantify the strength of the first order transition near to the 
critical Higgs mass several quantities have been considered. Measuring the dis- 
continuity of the Higgs condensate is practically tantamount (in three dimen- 
sions) to calculating the latent heat. Both quantities are easy to obtain but are 
very sensitive to finite volume effects. It is even more difficult to measure the 
interface tension for weak transitions. This we have observed in our recent work 
H where we have estimated the interface tension for a Higgs mass near to its 
critical value. 

In the literature mainly three methods are discussed to extract the interface 
tension from Monte Carlo studies. In the present work we shall employ and 
critically discuss the tunneling correlation length method flq, J7J , which has been 
used so far only in the analysis of 4-dimensional simulations H M in the context 
of the electroweak phase transition. 

In section 2 we define the model and discuss the methods used to extract the 
interface tension. The results obtained by the tunneling correlation method are 
presented and critically discussed in section 3. They are put into perspective 
with all available lattice measurements for arbitrary Higgs masses in section 4, 
which contains also our conclusion on the capability of the method. 



The Model and How to Measure the Interface 
Tension 



The model under study is a 3-dimensional 5'C/(2)-Higgs system with one com- 
plex Higgs doublet of variable modulus. The gauge fields are represented by 
unitary 2x2 link matrices U x>a and the Higgs fields are written as & x — p x V x . 
Px = iM^x ®x) i s the Higgs modulus squared, V x an element of the group 
SU{2), U p denotes the SU(2) plaquette matrix. The lattice action is 



S = /? G ^(l-Itrf/ p )-/3 ff Y, ^r($+C/ XiQ $ x+Q ) 

P ' l=(x,a) 

+£G£ + A*G£-i) a ) a) 



(summed over plaquettes p, links I and sites x), with the gauge coupling (ic, 
the lattice Higgs self-coupling /3r and the hopping parameter /?#. To relate the 
lattice couplings to continuum quantities the notations of la] are used. 

Let us recall the procedures used for extracting the interface tension of the 
electroweak phase transition. The first one, the two-coupling method, requires 



a careful extrapolation to the limits (3hc ±£ of simulations in two coupled 
subvolumes ||-Q supposed to be kept in the two phases. The variants differ 
in the way this limit is taken. In our version of the method O] we have used 
a simultaneous multihistogram technique for a system consisting of two parts 
in order to estimate the free energy difference between the homogeneous and 
mixed states at Phc- 

The majority of results for the interface tension has been obtained with 
the second method based on the two-state signal in the histogram of an or- 
der parameter like quantity o [[10], [13], ||, [|. Here the whole system is simu- 
lated at the pseudocritical point. The bulk variable under consideration is, say, 
the volume average of the modulus squared o = -^Ul{§ + <&) or the average link 

o = itr($+t/$). Usually, the interface tension a is estimated comparing the 
minimum and the two maxima of the doubly peaked histogram p(o) by 
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where L is a typical linear extension of a surface eventually dividing the system 
into different phases (the smaller extensions of a cylindrical system), a denotes 
the lattice spacing and T c the corresponding transition temperature. Different 
entropy factors in the thermodynamic weight of the mixed state have been ig- 
nored for simplicity in (|2|), but are necessary to extract a reasonable estimate 
for the interface tension, in particular if histogram data of lattices with various 
geometries are used simultaneously |1J] . Applying a formula like (g) one tacitly 
assumes that (i) the order parameter o is well-chosen in order to yield a his- 
togram with clearly separated maxima and a broad minimum in between, (ii) 
the case of equal height of the maxima is near to the phase equilibrium, (in) 
a minimal surface spanning through the lattice separates the pure phases from 
each other. 

Lattice studies of the electroweak transition in the case of the Standard 
Model for realistic values of the Higgs mass are known to possess a transition 
which is very asymmetric and weak. Therefore the conditions (i) to (Hi) are 
hardly fulfilled. By "asymmetric" we mean that the fluctuations of suitable 
order parameters o are stronger in one (the Higgs) phase than in the other 
(symmetric) phase. This fact makes it more and more difficult to use histogram 
methods based on (g) when approaching the critical Higgs mass. 

In order to cope with overlapping histograms we have recently proposed a 
method to split histograms p(o) into pure phase and mixed phase contributions 
[pi even under realistic circumstances of a very weak transition. This has enabled 
us to employ the histogram reweighting technique to find the pseudocritical 
coupling by the equal weight criterion. As a by-product, we have also obtained 



the latent heat from the variation of both pure phase thermodynamic weights 
near to the transition point and the thermodynamic weight of mixed states. 
Extracting, however, the interface tension from the relative weights of pure and 
mixed phases still depends on the simplifying assumption (Hi) above. Collecting 
data from lattices with different aspect ratios and extrapolating in the smallest 
linear extension to infinity we |5| obtained at the physical Higgs mass rriff = 
64.77 GeV (which corresponds to T c = 150.9 GcV) for the S?7(2)-Higgs theory 
without fermions an estimate of a/TJ? = 2.1 x 10~ 4 . 

All methods discussed so far are focusing on changes of (volume or sub- 
volume) averaged variables o in mixed phase systems with minimal interfaces. 
There is a third method p, uj to deduce the interface tension a from a tunnel- 
ing correlation length £tunnci a t the phase transition. To be more precise, it is 
the dependence of this correlation length on the geometry of the system which 
allows to extract a. A first test of the credibility of the method in the case of 
the electroweak phase transition has been undertaken in Ref. far from the 
critical Higgs mass and in Ref. M at even smaller Higgs mass, both within the 
4-dimensional framework. 

In contrast to measurements of the temperature dependent Higgs mass on 
both sides of the transition J5J (which requires to separate pure phase samples) 
the emphasis is here on correlations due to different phases in coexistence. The 
tunneling correlation length is measured in a very elongated volume L 2 x L z , 
stretched along the z-direction. In 4-dimensional simulations the additional 
Euclidean temporal extent is understood to represent the (inverse) temperature 
and cannot be modified in practice. The connected correlator C conn (zi — z<i) of 
o(z) = ^2 X S X3Z tr ($x^x) between two equal- z slices defines the correlation 
length ^tunnel which is expected to vary with the transversal extent of the lattice 
as 



^tunnel 



oc exp (a 3 A) , A = (aL) 2 . (3) 



The interface tension of the original 4-dimensional theory a is related to a 3 
through a — a 3 T c . 

Using a semiclassical expansion for a scalar Lp A theory in id including qua- 
dratic fluctuations around a kink solution J7J , the mass (energy) gap m = ^tunnel 
has been calculated (in lattice units) as 




(4) 
with no additional L dependence in front of the exponent. Therefore, expressed 



in terms of the dimensionless parameter 
aL 2 a 2 



Tr 



(5) 



the following scaling behaviour is expected to hold for the tunneling correlation 
lengths (if they are measured in units given by the transverse lattice size) 

mi a t L = Cy/x exp(-x) . (6) 

Strictly speaking, this result is valid for systems within the universality class 
of the Ising model, but it has been confronted with the 4-dimensional SU(2)- 
Higgs model in Ref . M . In this first application in the electroweak context (at 
smaller Higgs mass and, consequently, higher a) it has been shown that the 
simple perturbative one-loop result (0) is reached from above with increasing 
transverse lattice size and is valid only beyond x ~ 1 (cf. Fig. below). It 
has been argued already in [lq ] that higher order corrections may lead to a 
pre-exponential power in L with an exponent different from zero in M) R. Con- 
centrating on the roughening of the interface in a capillary wave model beyond 
the Gaussian approximation [17J] it has been found that (Q) gets a correction 
factor (l + 0(T c /(aL 2 a 2 ))). 

Recall also that in all derivations it has been assumed that the tunneling 
correlation length is much larger than the typical correlation lengths in the pure 
phases. In our recent studies at a physical Higgs mass of roughly 65 GeV || 
we have measured the Higgs correlation lengths near to the phase transition in 
the pure phases as £brokon/a = 13.40(41) and ( symm /a = 9.71(29) which are not 
so small compared to the tunneling correlation length l/(mi a to) as one will see 
later. Furthermore we have visualised there a typical mixed-phase configuration 
which had very rough interfaces separating different phases. 

Therefore, in order to extract an interface tension at this very weak first 
order transition (with very rough interfaces) we assume that the tunneling mass 
gap can be parametrised in a more general form as function of the transverse 
extent L 

m^t i = ci 7 exp (- ai a t L 2 ) , (7) 

with the fit parameters c, 7 and ai at . Using the 3d continuum gauge coupling 
g 2 = 4/(/3(ja), we can put a\ at into relation to the 3d dimensionless interface 

1 This has been confirmed in a recent two— loop calculation Jl™. 
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As a check of this assumption we have to compare the interface tension with 
the result of other methods when these are available. 



3 Results and Discussion 

The lattice model (||) is used as in ||, in particular with the same update 
algorithm as described there. We are dealing with the two cases corresponding to 
the Higgs masses M# = 70 and M* H = 57.4423 GeV (denoted in the following by 
57 GeV). These cases correspond to X s /gl w 0.095703 and A 3 /flf| « 0.0644457, 
respectively, (A3 is the 3d continuum Higgs self coupling) via 



- = -(^M 2 0) 
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The Higgs mass Mjj (in GeV) differs numerically only slightly from the physical 
Higgs mass m# in the Ad theory without top. 

The smaller Higgs mass is chosen in accordance with the work [g_3| (in their 
notation referred to as m* H = 60 GeV). Correlation function measurements 
have beeN Taken after each 10th iteration. The maximum of the integrated 
autocorrelation time for this quantity was about 26 at the smaller Higgs mass 
and 16 at the larger Higgs mass (each in case of the largest measured transverse 
size). 

To obtain the tunneling correlation length ^tunnel requires some tuning. At 
first, the appropriate hopping parameter value j3u has to be tuned, separately 
for each transverse size of the system while keeping near to the bulk critical 
value Phc, to the maximum of the tunneling correlation length. The actual 
longitudinal size of the lattice has been chosen three to four times larger than 
the correlation length one is going to measure. 

In the tables we quote the statistics for all lattice geometries L 2 x L z at the 
respective (5h (having the maximal tunneling correlation length) and its corre- 
sponding value m\ a t. We indicate also the inverse transverse correlation length 
wj_ iat and the total statistics which went into the search for the minimum. 
Additionally, we have checked at the larger Higgs mass for geometries 12 2 x 96 
and 20 2 x 128 that mi a t does not change within the errors for larger L z . 



L 2 x L z 


TOlat 


# msmts 


0H 


TOXlat 


total # msmts 


4 2 x 32 


.2180(106) 


8000 


.343600 


— 


36000 


6 2 x 64 


.1543(117) 


10000 


.343000 





24000 


8 2 x 64 


.1165(081) 


6000 


.342700 





8000 


10 2 x 64 


.09150(531) 


6000 


.342700 


.1008(29) 


10000 


12 2 x 64 


.07060(339) 


6000 


.342694 


.0782(21) 


12000 


14 2 x 96 


.05200(380) 


10000 


.342700 


.0632(21) 


19000 


16 2 x 128 


.04195(286) 


10000 


.342688 


.0509(15) 


21000 


18 2 x 128 


.03517(331) 


12000 


.342686 


.0425(14) 


20000 



Table 1: 



Statistics at M* H = 57 GeV 



L' 1 x L z 


raiat 


# msmts 


I3h 


m±i&t 


total # msmts 


4 2 x 128 


.276(6) 


3000 


.345000 


— 


18400 


6 2 x 64 


.1933(044 


10000 


.344200 


— 


60000 


8 2 x 64 


.1489(053) 


20000 


.343800 


— 


40000 


10 2 x 64 


.1174(043) 


25000 


.343540 


.1345(17) 


80000 


12 2 x 64 


.09530(278) 


50000 


.343540 


.1080(13) 


110000 


14 2 x 64 


.08031(267) 


25000 


.343540 


.0889(16) 


25000 


16 2 x 64 


.06910(154) 


25000 


.343540 


.0761(09) 


40000 


20 2 x 64 


.05367(163) 


15000 


.343560 


.0599(12) 


51000 



Table 2: Statistics at M 



H 



70 GeV 



In order to extract the correlation length funnel we first check by inspection 
whether the local mass has a plateau. We define a local lattice mass m\ a t(z) 
at correlation distance z through a fit of three subsequent values C conn (z — 1), 
Cconn^) andC conn (z+l) of the correlation functions of o(z) = J2 X S X3 ^ Z tr ($+$3.) 
to a hyperbolic cosine shape A (exp(— mi a t(z) z) + exp(— m\ &t (z) (L z — z))). This 
ansatz is essential in order to observe a plateau in the local mass versus z. In 
Fig. |l| 77ii at (z) is shown for one particular example. We contrast this with a 
local mass defined through a single exponential fit to the three neighbouring 
values of C CO nn(z) above. No plateau at all can be identified using the latter 
definition of a local mass. Finally, the (inverse of the) correlation length ^tunnel 
is obtained by a global hyperbolic cosine fit over the plateau range that we have 
identified. The behaviour of the inverse correlation length near to its lowest 
value is demonstrated in Fig. 0. 

In Fig. |3J we show the inverse tunneling correlation length multiplied by the 
transverse extension, mi at L as function of the transverse lattice area L 2 for both 
values of the Higgs mass under discussion. The strongly different exponential 
slopes reflect the weakening of the phase transition. The curves correspond to 
a least square fit with the ansatz (Q). 
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Figure 1: Example of the local lattice mass at Mjj = 57 GeV and 16 2 x 128 as 
function of z 
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Figure 2: Inverse correlation length at MJj = 70 GeV and 12 2 x 64 as function 
of(i H 



From this fit we extract the dimensionless 3d interface tensions 03 j g\ = 
0.0224(56) for the lower Higgs mass M* H = 57 GeV, and a 3 /g£ = 0.0049(18) 




500 



Figure 3: Fits for the inverse tunneling correlation lengths at Mjj = 57 GeV 
(circles) and M^ = 70 GeV (squares) 



for M%j — 70 GeV. The first number is in very good agreement with the value 



presented in Fig. 15 of 



|13j (and - with the use of eqs. 
= 0.0217(22). 



(11.6) and (2.8) there 



translated into) ct^/g^ 

Using these numbers we find the interface tensions a 



— = 3.24(80) x 10-\ 



M* H = 57GeV 



(10) 



and 



o 



^3 = 70(26) x 10- 5 , 



M* H = 70GeV 



(11) 



corresponding to the 4-dimensional model without top quark. Taking into ac- 
count the top quark the numbers change to 3.57(89) x 10~ 3 and 77(29) x 10~ 5 . 
The relations between 3-dimensional and 4-dimensional quantities have been 
obtained as in || using (in the same notations as there) table || which is recalcu- 
lated here for M^ = 57 GeV. Though the top contribution apparently changes 
the interface tension only insignificantly (A g is small) some of the fermionic 
one-loop corrections to the 4d couplings are already too large and the other 





M* H = 57 GeV 


M* H = 57 GeV 




without fermions 


with top 


A 9 


-0.01322 


6.712 x 10~ 5 


A A 


-0.02383 


1.494 


K 


-0.02938 


0.9443 


Ay 




-0.01301 



Table 3: Corrections in eq. (6) of [5] 



physical numbers should be taken with great care. The reported pseudocriti- 
cal (3h values for the largest transverse sizes from tables [l] and g are related 
to critical temperatures T c — 127.2 and 150.1 GeV and physical Higgs masses 
rriH = 52.35 and 64.55 GeV of the 4d theory without top for the lower and 
larger M^. The four dimensional MS gauge coupling g 2 (mw) has the value 
0.423 and 0.422, respectively which is close to that of the standard model. 

Our present result for the case of M^ — 70 GeV is larger by a factor 3.3 
than the previous rough estimate |j|. The latter was obtained as a result of 
our equal weight histogram method, finally based on a global infinite volume 
extrapolation of the mixed phase thermodynamical weight for lattices of various 
aspect ratios. 

In the present fits we obtain effective exponents 7 = 0.25(11) and 0.095(45) of 
L, respectively, which become smaller with decreasing strength of the transition. 
Fixing the exponent to 7 = 1 as suggested by eq. (0) we would be able to present 
only a local fit to the few highest transverse areas L 2 . The interface tensions 
would be estimated by this fit as follows: as/gf = 0.023 for Mjj = 57 GeV 
and 0.016 for M^ = 70 GeV. In the first case this would be still acceptable 
comparing with the result of the Helsinki group [H3] for that Higgs mass. But 
the interface tension evaluated at Mjj = 70 GeV in this way does not follow the 
general trend of the interface tensions from d = 3 simulations which have been 
collected in Ref. Q. 

In view of the arguments above, we consider the fit with the free ansatz 
(0) more serious than the results of the fit confined to 7 = 1. We recall that 
also in Ref. || the latter fit has been successful only at large values of the 
scaling variable x — aL 2 a 2 /T c . The one-loop scaling law (g|) was meant to 
hold independent of the particular system, irrespective of the actual value of 
the interface tension. With the values of a obtained now and in Refs. [g, [| the 
underlying mass gaps for various Higgs masses and transverse lattice sizes seem 
to follow another universal law if expressed through the scaling variable x 



miat L = C' exp(-x) 



(12) 
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Figure 4: Scaling law for the mass gap followed by data at ttih — 35 GeV (open 
squares) [9], ttih = 49 GeV (open circles) [8], Mjj = 57 GeV (full circles) and 
M* H = 70 GeV (full squares) [this work] 

instead of (|lj), except for the smallest transverse extensions L in each case. 
This can be seen in Fig. [|. In this figure we also show the asymptotic behaviour 
(|) with C = 1.352 as predicted by @, [uj. 

There is one warning in the data concerning the applicability of the analytical 
formulae used throughout the literature in order to extract the interface tension 
from lattice data of the tunneling correlation length. As mentioned above it is 
implicitly assumed that the correlation lengths of the two phases in equilibrium 
are small compared to the tunneling correlation length. From this point of view 
the present data on m\ B ± for Mjj = 70 GeV (where we have measured Higgs 
masses separately for both phases at the transition temperature) indicate that 
transverse sizes L < 10 would have to be excluded from the analysis. 



The other concern is caused by the fact that the transverse correlation length 
£j_ does not decouple from the transverse lattice extension at small L and keeps 
growing for all transverse sizes considered. For instance, at M^ = 70 GeV, the 
transverse correlation length amounts from 75 to 83 per cent of L on our lattices 
with transverse sizes L > 10 which are, on the other hand, the only acceptable 
ones in view of the criticism discussed before. In the case of M^ = 57 GeV the 
transverse correlation length is even larger compared with the transverse size L 
(100 to 130 per cent) on lattices with L > 10. It should be mentioned that the 
ratio of the transverse correlation length to the transverse size is consistent with 
the ratio of the bulk correlation length to the system size for cubic symmetries 



11 



on top of the phase transition (measured without separating the Monte Carlo 
sample into pure phase configurations). 



4 Overview and Conclusions 



To compare our 3d results for the interface tension with those of Ad measure- 
ments we follow the procedure outlined in pQ|. The measurements in the Ad 
approach have been performed at a different gauge coupling. The measured 
renormalised gauge couplings do not seem to change significantly with the Higgs 
mass in the so far reported region from 18 to 49 GeV M, ElJ and vary from 
0.56 to 0.59. For simplicity (and due to missing calculation) it is assumed as 
in [ pp| t hat the measured renormalised gauge coupling roughly corresponds to 
the MS running coupling. It is then different from the value g 2 {mw) = 0.42 
corresponding to our calculations described in the last section. 

Using the prescription to relate 3d and 4<i parameters || we calculate Higgs 
masses, critical temperatures and interface tensions for increasing Ad gauge cou- 
plings keeping the respective X^/g 2 fixed. The numbers are collected in table [|. 
Note that the Higgs mass is slightly moving to lower values while the ratio a/T^ 



M* H /GcV 


g 2 {m w ) 


mn/GeV 


T c /GeV 


«/T c 3 


57 


0.423 


52.35 


127.2 


0.00324(081) 




0.560 


50.86 


108.2 


0.00539(134) 




0.570 


50.75 


107.1 


0.00556(138) 




0.580 


50.64 


106.0 


0.00574(143) 




0.590 


50.54 


104.9 


0.00592(147) 


70 


0.422 


64.55 


150.1 


0.00070(26) 




0.560 


62.99 


128.1 


0.00115(43) 




0.570 


62.88 


126.8 


0.00119(45) 




0.580 


62.76 


125.5 


0.00122(46) 




0.590 


62.65 


124.3 


0.00126(47) 



Table 4: Higgs masses, critical temperatures and interface tensions for various 
Ad running gauge couplings 



becomes much bigger, largely due to the critical temperature getting smaller. 

Having the caveats of the last section in mind we compare now in Fig. @ our 
interface tensions a with those measured by different methods in the Ad theory at 
various Higgs masses as function of the physical Higgs mass mn ■ In order to do 
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Figure 5: Collected results on the interface tension as function of the Higgs mass 



this we read for our 3d data (full points) the corresponding Higgs mass (without 
fermion contributions) from table 0, for example at g 2 (mw) = 0.58. This value 
has been deduced as mentioned above from the measured renormalised gauge 
coupling. 



In the 4<i data the Higgs masses are either measured [jT(J, || O, |[ 12 (no 
errors are taken into account in the horizontal scale) or estimated |j7 The agree- 
ment between the various approaches in 3d and 4d is remarkable if the 3d and 
Ad parameters are correctly mapped onto each other. This nicely demonstrates 
the validity of dimensional reduction in the considered Higgs mass range. 

Our new data points are the following ones, expressed in terms of the 3- 
dimensional and the 4-dimensional interface tensions. We find at Xs/gl = 
0.0644457, g 2 {m w ) = 0.423 and (3 G = 12 {M* H = 57 GeV) 



2J = 0.0224(56), ^ = 3.24(80) x 10~ 3 

and at \ z jg\ = 0.095703, g 2 (m w ) = 0.422 and /3 G = 12 {M* H = 70 GeV) 
a 3 



(13) 



- 0.0049(18), — r = 70(26) x 10 

93 T c 



-5 



(14) 
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We observe an approximate scaling law expressing the energy gap for all 
Higgs masses in terms of the dimensionless variable x (see eq. (||)) without the 
prefactor ^/x which had been suggested by one-loop perturbation theory for 
interfaces in the case of binary systems. 

We emphasise that the tunneling correlation length method works well even 
near to the critical Higgs mass where other methods relying on discrimination 
of histogram peaks and minima are already difficult to apply. 
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